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Attention: 

Subject: 


AP13-F; Ray Weems;' Contract Manager 

Contract NAS8- 26642, "Study of Multilayered Insulation 
Pipe Penetration," and "Thermal Acoustic Oscillation," 
Quarterly Progress Report 


Report Period: 26 July through 25 October 1974 
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This letter is submitted in accordance with the provisions of Exhibit C of the Y* 
subject contract and describes the direction of effort and accomplishments for^» 
the period shown above. This study is being performed by personnel in the ^ 
Thermodynamics & Structures Section of the Lockheed-Hunts ville Research & 
Engineering Center. -The Contracting Officer^ Representative (COR) for this ^ 
contract is Mr. E. Haschal Hyde, EP43. r 
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Study of Multilayered Insulation Pipe Penetration 
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The agreement between the analytical predictions of the net heat leak to a 
stored cryogen (L^) caused by 4 a metal penetration through the blanket of 
multilayer insulation and the experimental results as modeled was shown in 
progress report LMSC-HREC PR D390395 to be excellent. For brevity, a 
single example is shown again in this report for discussion purposes. ..The , 
~~ example shown (Fig. 1) is for a 7.62 cm o.d. stainless steel pipe with a 'wall 
thickness of 0.216 cm. The penetration was insulated over a 91-4 cm length 
_ _ with a 2.54 cm thickness blanket. The cryogen storage' vessel is“aluminum~ 

with the exception of a stainless steel access port on the opposite side from 
the penetration. All penetrations tested were stainless steel except for two 
fiberglass overwrapped Inconel units. The penetrations are mounted on 
individual stainless steel flanges which are bolted to the cryogen tank. 
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Figure 1 shows that the mating surface temperature equilibrated at approxi- 
mately 86°K instead of the temperature of 77.4°K (the cryogen tank is 
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maintained at a constant 1 Oi atmosphere) The differential in the equilibrium 
temperature is caused by a contact resistance between the flange/cryostat mount- 
ing surface; Extensive efforts were made to minimize this contact resistance, 
with the Results being in the range of 5 to 9°K. This parameter has a direct 
effect on the total efficiency of, the system and the resultant boil-off rate of the 
cryogen. ’ , 

Other parameters which must be considered in the analysis of the penetration 
systems are the end temperature of the penetration, the diameter and wall thick- 
ness of the penetrations, the material used for the penetration, the insulation 
length and thickness and the cryogen being stored. 

The results shown in Fig. 1 indicate very well the capability of the analytical 
model as how programmed to predict the system being .used as the experi- 
mental verification of the model. Due to the limited flexibility as planned for 
the computer program for checkout purposes, variables including cryogen and 
penetration material properties were not originally coded into the program. 

That effort is now being undertaken with initiation of the final results for all' 
these variables anticipated to begin during the November-December perform- 
ance period. The computer program will then be executed for the following, 
variables. \ „ - - * 

' v * * ; - 

• Cryogen* 

• * Penetration material 

• Penetration insulated length 

, • Penetration insulation thickness 

• Contact resistance 

« ' • Penetration temperature at end of penetration 

insulation ’ 

• Penetration diameter . * ♦ 

• Penetration wall thickness 1 


• Penetration radiation properties 



The cryogens to be investigated include* helium, hydrogen and oxygen. The 
cryogen storage, container will at all_times bejasj3umed to.be aluminum. „The_ 
penetration materials analyzed will include stainless steel, fiberglass and 
Inconel. Insulation lengths and thickness will be varied over typical applica- 
tion ranges, as will the diameter and wall thicknesses of the penetrations 
analyzed. The contact resistance will be varied between the limits of 0°K 
drop (limiting case) to 20°K. The temperature of the penetration surface at 
the end of the insulation blanket will be varied depending on the insulation 
length and cryogen, with a judgment value' being used which will be based on 
the results of the analytical/experimental verification program and the cryogen. 
This is necessary because of the large range in stored cryogen temperatures 
being considered (helium = 4°K, hydrogen = 20?K and oxygen = 90°K). 
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The surface finishes ol the penetrations will be assumed to be polished with 
the exception of the fiberglass. Typical radiation properties of these materials 
for polished (metal) or "as produced” (fiberglass) will be used. Temperature- 
dependent conductivity values will also be incorporated. 

* • . j 

The computer program is presently being modified to include the necessary 
properties of the various materials in either equation or table look-up form, 
depending on that property required. 

Results of these analytical investigations, which will ultimately lead to a hand- 
book form of penetration results, will be reported on as they become available. 

• Thermal Acoustic Oscillations Analytical Study / 

The purpose of this study is to develop analytical models for predicting the 
effect and magnitudes of thermoacoustic* oscillations" on the penetration heat 
leak to cryogens and ^experimentally verify these models. The oscillations 
develop as a- result of large thermaLgradients imposed on a compressible fluid. 
The expansions and compressions of the fluid generate pressure waves in much 
the same manner as pushing a piston through a gas-filled pipe. These thermo- 
acoustic oscillations can greatly enhance the heat losses over those predicted 
by pure conduction theory. The Lockheed-Huntsville Thermoacoustic Oscilla- 
tions (TAO) computer program has been modified for use in this study. In , 
addition, a computer program for performing a power spectral density (PSD) 
analysis is being utilized. ' r f, : *f t j, “ 

Progress during this quarterly period is summarized in^the following categories: 
(1) statistical analysis; (2) TAO parametric study; and (3) power spectral density 
calculations. > , . "• 

i ” t 1 

Statistical Analysis- The primary objectives of this statistical study is to 
determine the predominant amplitudes and frequencies of the thermal acoustic - 
^oscillations. The analysis involves first the determination of statistical' pa- 
rameters such as mean value, variance and standard deviation of the dependent 
variables (i.e. , pressure, temperature,, density and velocity profiles). In addi- v 
tion, several important statistical functions are determined and analyzed. These 
-include the auto -correlation function, the power spectral density function, the 
cross -correlation, function, and the cross spectral density function. The follow- 
ing' is a brief discussion of the definitions and significances of the statistical 
parameters and functions (Ref. 1)\ t w 

^ *• S , ** * ^ S. I 

•We'start with a time series x(t) having oscillatory properties and being valid 
a sufficiently long period 0<t<T.. The series x(t) may represent the cal- 
culated pressure, temperature, density or velocity solutions of the thermal 
acoustic oscillations computer program.' .The statistical parameters of x(t), 
namely, the mean value x, the variance <y the mean square value > 2, anc j 

the standard deviation a , are, defined as follows: 
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x = T J 

f x(t) dt 


(1) 

O 

/* T 2 



"x --Tj 

f [x(t) - x] dt 

. 

(2) 

C 

f T 



< -+J 

f X 2 (t) dt 


(3) 

o 




** -s°l 


(4) 


The mean value is simply the average of all values within the time interval 
of 0 < t < T. The variance is the mean square value about the mean. The 
variance is a dynamic or fluctuating component of the time series. It is a 
measure of the degree of fluctuation. The mean square value represents the 
general intensity of the time series. It is the average of the squared value of 
the time history. The standard deviation is simply defined as the positive root 
of the variance. , 

By expanding Eq. (2), it can be shown that the variance is equal to the mean * 
square value minus the square of the mean value. That is, 

"i = +* - < s > 


Auto -Cor relation Function : The auto-correlation function of an oscillatory 
time series x(t) describes the general dependence of the values of -x at one 
time or the value at another time. Consider a time history record x(t) illus- 
trated in Fig. 2. An estimate for the auto -cor relation between the values of 
x(t) at times t and t + x may be obtained by taking the product of the two values 
and averaging over the observation time t. The observation time , t is supposed - 
to be sufficiently large to co'v’er" many times the periods of the waves Contained 
in x'(t). In equation form, the auto-correlation function, R (r), is defined as 

- -- - 1 

T 

x(t)x(t+T)dt (6) 

The quantity R x (r) is always a real-valued even function with a maximum at 
T = 0, and may be either positive or negative. In terms of the auto-correlation 
function, the mean value of x(t) is given by 


R x < T > 


= T / 
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In words, the mean value of x(t) is equal to the positive square root of the auto-' 
correlation as the time displacement becomes very long. Similarly, the mean 
square value of x(t) is given by 

+ x 2 =R X (0) (8) 

That is, the mean square value is equal to the auto-correlation at zero time 
displacement. " • 

Power Spectral Density Function: The power spectral density function (also 

called the auto- spectral density function) of x(t) describes the general fre- 
quency composition of x(t) in terms of the spectral density of its mean square 
value. The" mean square values of the time series x(t) in a frequency range 
between’f and f + Af, i.e., x(t, f, Af), may be obtained by filtering the series 
with a band-pass filter having sharp cutoff characteristics, and computing the 
average of the squared output from the, filter. This average squared value will 
approach an exact mean square value as the observation time, becomes large. 
That is, for sufficiently large t, 1 - 


(9) 

Where \\j ^ (f, Af) is the mean square value of x in the frequency range from 
f to f + Af. 

The power spectral density function of x(t) is defined as 


* (f,Af) 


■+P 


(t, f, Af) dt 


G (f) 
x' 


lim < f ' 

Af-*- 0 • Af 

£ ^ t f * 2(tf: f - Af )] dt ” 


^ “( 10 ) 


Hence, the power spectral density of x(t) at the frequency f, i.e., G (f), is the 
mean square value of x(t) at f per unit frequency. 

An important property of the power spectral density function lies in its relation- 
ship to the auto-correlation function, R x (t ). For a set of stationary oscillatory 
data, the two set statistical functions G (f) and R (t) are related by a Fourier 
transform as follows 
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G (f) 
x 


Where j = 

Other statistical properties of x(t), such as mean value and square value, are 

also related to G (f)* , % - ‘ 

x , , - , « *.* 

A typical plot of power spectral density function versus frequency is shown in 
Fig/ 3. The figure shows that the original time series x(t) contains two pre- 
dominant waves of frequencies fj and f^* - / 

The principal application for a powe r" spectral density function measurement 
of physical data is to establish the frequency composition of the data which, 
in turn, bears important relationships to the basic characteristics of the phys- 
ical system involved. For our present study of thermal acoustic oscillations, 
a power spectral density analysis of the analytical results of t pressure, tem- 
perature, etc. , will enable us-,to determine the predominant frequency or fre- 
quencies of the system. ' - v * ; 



Fig. 3 - Power Spectral Density Plot 


= 2 f R x (t) e " j2 n { * 

- 00 

= 41 R x (t) cos 27T ft dt 


( 11 ) 
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TAP Parametric Study : Computer runs which are being made consist of two 

types: (1) fixed length-to -diameter (L/D) tubes with parametric acoustic 
Reynolds number for use in constructing an existence diagram, and (2) gen- 
erating data tapes with complete flow fields versus time for use in the PSD 
program. Two curves of particular interest are presented in Figs. 4 and 5. 
These are plots of temperature and pressure versus time for an open tube 
with L./D = 300. These are shown to illustrate the significance and need for 
a PSD analysis. Note that both profiles show a sequence of "high-frequency" 
waves superimposed on a "lower frequency" wave. The PSD analysis will 
operate on this type of data to extract the dominant frequencies of the waves. 
The temperature plot is shown at z = 0.5 (near warm end) and the pressure 
plot is shown at z = 0.95 (near cold end) which illustrates that this multi-wave 
characteristic is present at both ends of the tube. 

One of the primary objectives of this study is to determine conditions which 
produce thermal acoustic oscillations. This result is being sought by per- 
forming a parametric analysis using the TAO program. 

The parameters being studied are those given by Von Hoffmann et al. (Ref. 2): 
(1) ratio of hot to cold temperature (T^/T^); (2) length -to-diameter ratio of 
tube, L/D; (3), acoustic Reynolds number. Re; and (4) ratio of total tube length 
to cold part length, L/L . Figure 6 is a reproduction of the Von Hoffmann 
figure with TAO cases shown for comparison. Several other TAO cases are 
being run in order to construct a curve similar to that of Von Hoffmann. The 
data points shown are in reasonable agreement but with several apparent dis- 
crepancies. For example, at Th/,T C = 10, the TAO solution shows that oscil- 
lations occur for values of the dimensionless parameter A = 100 while Von 
Hoffman shows no oscillations above A = 30. This is most likely due to the 
nonlinear effects included in the TAO solution while Von Hoffman's analysis is 
based on linearized equations. 

The solid Line indicates the trend of the curve for defining the existence 
criteria for thermal acoustic oscillations. Several other cases are being 
processed to refine this curve. This portion of the study will be completed 
during the November -December period. 

Power Spectral Density Calculations: The Lockheed-Huntsville Statistical 

Analysis Computer Program is being utilized in this study. This program 
is~a multi -function generation program for~statisticaTanalysis , ~i.e it * ~~ ~ 
computes more than just the spectral density function. 

The inputs x, y can be any functions of time such as temperature and pressure 
These are supplied via magnetic tape from the TAO solutions. The program 
then calculates mean values, mean squared values, standard deviations and 
variances. The data are then detrended and tapered to give zero mean value 
and slope. The Fourier transform of the complex function z = x + jy is now 
computed using a Fast Fourier Transform (FFT) algorithm. The first calcu- 
lations are then made for a number of statistical functions. The PSD function 
is of primary interest here but the cross correlations may prove useful also.* 



Temperature, T/T 
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Fig. 4 - Temperature Profile at 
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The program is operational for the general cases of inputting the x(t) via cards. 
Due to the complexities and large amounts of TAO data, the program has been 
modified to read the TAO solutions via magnetic tape. 

During this performance period, the PSD program was run for several cases. 
The objective of these test runs was to verify the operation of the code modifi- 
cations. The first test case consists of two square wave signals as shown in 
Fig. 7. The power spectral density function corresponding to the input signals 
is shown in Fig. 8. This basic checkout proved successful. 

The next checkout case consists of a portion of temperature and pressure 
wave calculations from the TAO program'. The input signals are given in 
Fig. 9. The top curve is temperature and the lower one is pressure (dimen- 
sionless). Figure 10 shows the calculated power spectral density. The pres- 
sure input here is a portion of a pure sine wave and thus the power spectral 
density shows a "spike' 1 at the one frequency and zero elsewhere. This is 
correct for a pure sine wave. The temperature PSD plot shows the content 
decreasing at high frequencies. This is again f he proper behavior. ’ 

A third case has been processed which consists of a larger portion-of an actual - 
TAO solution. Figure 11 shows the temperature (top) and pressure (bottom) 
input signals. Figure 12 is the corresponding power spectral density plot. 

This case appears to be correct, but further analysis is being done to verify 
the proper operation for this case. When this is finished, the TAO tapes for 
various L/D ratios will be run. These will be presented in the next progress 
report. 

• Thermal Acoustic Oscillations Experimental Verification Program 

( 

The experimental verification program for the thermal acoustic oscillations 
analytical predictions has been prepared for initiation. In order that a more 
permanent form of experimental data may be obtained than would have been 
possible using an osc lllosc ope/camera system, an FM recorder will be used 
to obtain a permanent record of the frequency/amplitude output of the pressure 
transducer. This will enable repeat playbacks of the data and will also allow 
an analog -to-digital data translation for input to the power spectral density 
prog ram. if time allows. _This portion of the program is desirable although not 
mandatory. The added capability will enhance the interpretation of the data if 
it is possible to input the experimental data into the PSD program. 

The initial quantity of liquid helium to conduct the tests has been ordered, and 
testing is scheduled to begin within two weeks. The 0.95 cm outside diameter 
tube has been selected as the "control" tube, and will be tested over a range 
of 50 < L/D < 300. Both stainless steel and aluminum tubes are to be tested. 
The complete test matrix of tubes and their characteristics are shown in Table 1. 

• Long Term Cryogenic Storage Techniques 

The basic system configurations for long term storage of liquid cryogen for a 
low-g environment have been selected. These configurations fall basically 
into four categories listed in Table 2. 
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Table .1 

TAO EXPERIMENTAL- VERIFICATION PROGRAM 
TEST TUBING SIZES 


— 

Tube i.d. 
(cm) ' 

Wall Thickness 
(cm) 

Tube o.d. 
(cm) 

l/d 

(-) 

L 

(cm) 

Stainless Steel Tubes. 

0.30 

0.089 

0.478 

150 

45 

0.658 

- . 0.147' 

- . Q.953 

50, 75, 100, 

33, 49.5, 66.1, 




150, 300 

98.8. 197.6 

1.297 • 

0.305 

1.905 

30, 150 

-38.8, 194.3 

Aluminum Tubes 

• 0.658 

0.147 

0.953 

50, 100, 150 

33, 66.1, 98.8 


Table 2 


SUMMARY OF CONFIGURATION CONCEPTS 

Configuration 1 

Single blanket of MLI 

Configuration 2 

MLI with vapor cooled shield 

Configuration 3 

MLI with solid cryogen shroud 

.3’ 

Solid hydrogen shroud 

.3" 

Solid nitrogen shroud 

.3’" 

Solid argon shroud 


Configuration 4 


MLI with radiator cooled fin 
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Fig 7 - Input Signal for PSD Checkout Run: Case 1 
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POWER SPECTRAL DENSITY 


JO0 NO AS 1008 , PAGE 2 
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Fig. 9 - TAO Sample Inputs for PSD Program: Case 2 
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Mathematical Models : .The analytical models for these configurations are 
based on on e -dimensional heat transfer through a homogeneous medium. 
There are two implications due to this basic simplification that are common 
to all configurations, namely: - 

• The cryogen tank wall is fully wetted 

• There exist no property or temperature gradients 
orthogonal to a radius. 


The latter implication impacts the modeling of such phenomena as penetration 
heat leaks, the heat transfer process in vapor cooled shields (Conf iguration 2), 
and the shorting of penetration heat loads at shields and cryogen shrouds (Con- 
figurations 2, 3 and 4) as will now be discussed. 


Penetration Heat Leaks ; The accompanying sketch (Fig. 13) indicates the method 
in which the penetration heat leaks to the liquid helium through support structures, 
vent tubes, electrical leads, and/or through the structure of the experiment itself 
will be modeled. The penetration load is the result of conduction heat transfer 
to the stored cryogen through several channels which are in parallel. In the ab- 
sence of heat shorts to vapor cooled shields or solid cryogen cooled shrouds, 
these heat penetrations through the parallel channels will be subjected to the 
same temperature difference, (T cr - ^ambl* In order to vary the penetration 
heat load in a consistent manner, a range of values will be assigned to a "pene- 
tration conductance." The conductance must characterize the heat penetrations 
through all penetration' channels and will therefore represent a summation of the, 
conductances through all channels. Thus, 



i = number of penetration 
channels (supports, 
vents, loads, etc.) 



Fig. 13 - Model for Penetration Heat Load 
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The penetration conductance is determined by:. (1) the geometry of the penetra- 
tion channel which may be characterized by A/L (where A may be taken to be 
dependent on L and an appropriate average value for A may be defined); and 
(2) the thermal conductivity of the penetration channel which may be taken 
to be temperature dependent to arrive at an average value for K. The 
heat load due to a penetration will be accounted for by parameterizing thermal 
conductance for the penetration heat load such that the range of conductances 
for various support and vent schemes are covered; A designer may then esti- 
mate the penetration conductance of his cryostat in order to select the appro- 
priate curve to be generated in this study. 

A preliminary estimate for the range of penetration conductances was based 
on using 347 stainless steel and fiberglass supports. The size of the, support 
members required to support cryostats weighing 45 kg and 227 kg for 10 g 
loadings with a safety factor of 2 were calculated. The lengths of the supports 
were taken £o be 0.3 m and 1.55 m along tKe penetration channel. Table 3 sum- 
marizes the influence of support material and liftoff cryostat mass on the pene- 
tration conductance. 


Table 3 ; ' . * . 

ESTIMATE OF THERMAL CONDUCTANCE FOR SUPPORTS 


, . > ^ f t 


- Material 

'* , 

' -4 ‘ 

Yield 

Strength 

V -2 

kg -m , 

347 Stainless 

9. Ox lO 7 

Steel C . „ 


Fiberglass 

9.1 x 10 6 

i 




L = 1.52 m 


0.47 x 10 


0. 146 x 10' 


L = 0.30 m 


1 . 1 72 x 10 


0.366x 10* 


The conductances listed in Table 3 do not account for penetrations through 
electrical leads and vent tubes. , One might double the conductances in Table 3 
to account for these effects. Also, a cryostat may be designed for retractable 
supports once the vehicle is placed in orbit in which case the conductances m 
Table 3 would \e larg er ; thsui The long term, conductance to the cryogen. To 
appropriately parameterize the conductances (see Table 5),. it appears that con- 
ductances should cover a range from 10“^ to 10"^ W/°R. 

Vapor Cooled Shield Efficiency : Vapor cooled shields immersed within the 

MLI are designed to allow the boiloff vapor to intercept a portion of the heat 
leak through the MLI and penetration channels. This implies the presence 
of a lateral temperature gradient within the shield which is in opposition . 
to the .one -dimensional heat transfer assumption on which the math model is 
based. To circumvent this difficulty an efficiency is defined forthe heat ab- 
sorbed by the boiloff vapor during its residence within the shield. Since the * 
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vapor cooled shield is essentially a heat exchanger, the efficiency of the 
shield is defined with the standard heat exchanger efficiency (see Fig. 14). 



C r yog en 


Fig. 14^ - Schematic of Vapor Cooled Shield 


The maximum heat is absorbed by the boiloff vapor if the vapor leaves at the 
temperature of the shield. In actual operation the vapor leaves at a tempera- 
ture below that'Of the shield. 
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Thermal.Short Efficiency: Efficient design of a cryostat often incorporates 

schemes for allowing vapor cooled shields and cryogen shrouds to absorb 
a portion of the penetration heat loads. This is referred to as thermal shorting 
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of the penetration heat load. The presence of a thermal short in a cryostat 
design violates the one-dimensional heat transfer assumption. To allow the 
inclusion of investigating the influence of thermal shorts in this analytical work, 
a definition for thermal short efficiency is made from which the magnitude of the 
thermal short may be determined. Reference should be made to Fig. 15. The 
maximum thermal short occurs if T p = T s h‘ In that case Qj may be defined for 
equal AK/L between shield/environment and shield/cryogen as 

'V* Sd.al = TT < T amb + T ery - 2 T S h> 



In actual^ operation Tp is above Tg^ ancLq^ is less than ^ideal* Thus, the 
thermal short efficiency is defined as the ratio of actual thermal short to the 
maximum attainable thermal short 


't.s. 


1 actual 

q ideal 


from . which T^ may be defined in terms of_ T^, 


cry’ 


T sh and 


as 


T 

P 


T , +-T 
amb cry, 


't.s 


■’(T* . + T • 2 T , : ) 

amb cry sh 
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The penetration heat leak to the cryogen is then 


q = - ~ (T - T ) 
• L cry p' 


and the magnitude of the thermal short is 


9 



amb 


+ T 
.cry 


- 2 v 


It should be noted that the magnitude of the thermal conductance is taken to be 
constant throughout the penetration channels in the above definition of thermal 
short efficiency. This assumption is probably not too restrictive because 
previous studies of cryogen storage systems, particularly shield designs, 
indicate that the optimum shield spacing is roughly equidistant between 
environment and cryogen. 


Configuration Equations : The equations used to describe a cylindrical cryogen 

.tank with spherical end caps was presented in the last progress report. Briefly, 
the heat leak through the MLI was given by the equation: 


q MLI = Zlr 

‘ 1 ■ 

where r , and r^, and T^ and T^ are the inside and outside cryostat radii and 
temperatures, respectively. is the conductivity of the MLI perpendicular 

to its grain, and L is the length of the cylindrical section of the container. 

The total heat load on the cryogen was given to be the sum of the heat leaks 
from-the experiment,- through the MLI, and through the penetrations i - The ^ * 

heat absorbing capacity of the cryogen was calculated from the equation for 
the cryogen mass. and its heat of vaporization which was assumed _to_be used 
uniformly over the storage time required of the cryostat. This heat absorb- 
ing capacity was then set equal to the total heat load and the equation solved 
for the cryogen mass. Finally a transcendental equation for r£ as a function 
or r j was found by calculating the mass of the cryogen based on its density 
and storage volume, and then eliminating the cryogen mass by equating this 
equation with the above described equation. Once ^ is found for a given fj 
the total system mass is calculated and stored. The value of r j is then incre- 
mented, the value of T2 is again determined and the resulting total system 
mass is compared with that previously calculated. This procedure is repeated 
until the minimum total system mass is calculated, at which time the minimum 
total system mass and the corresponding parameters are output. Figures 16 
through 19 show results obtained from the computer program using these 


2r l r 2 


(r 2 - r.) + in{r 2 /r . ) 


l 

/ 


k mli dT 


J T 



Minim urn Massr(kg) 
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Fig. 17 - Minimum Total System Mass for Hydrogen as a Function 
of the Length of the Cylindrical Section of the System for 
Con figuration 1 
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Hydrogen Storage 

Storage Time = 3 yrs 
Cylindrical Length = 3.05 m 
Ambient Temperature = 500°R 
Cryogen Temperature = 36.7°R 
Ambient Pressure =0 ^ 2 

Cryogen Pressure = 1.03 x 10 kg/m 
Experiment Heat Lpad = '0 ' * 

Safety Factor - 1 v ** 

MLI = Superfloe at 11.8 layers/cm 
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Tank + MLI 



Hydrogen Storage 

Storage Time = 3 yrs 
Cylindrical Length = 3.05 m 
Tank % Filled = 1 . . 
Ambient Temperature = 500°R 
Cryogen Temperature = 36.7°R 
Ambient Pressure ‘= O ' ‘ 2 

Cryogen Pressure = 1.03 x 10 kg/m 
Safety Factor = 1 

MLI = Superfloe at 11.8 layers/cm. 
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equations and procedure. The figures show the effect of storage time, length 
of cylindrical section, the percent which the cryogen tank is filled, and the 
experiment heat load on the minimum total system mass when hydrogen is 
stored using Configuration 1. The associated parameters are also shown on 
each figure. In particular, Fig. 16 shows the effect that the storage time re- 
quirements have on the minimum mass for, a spherical configuration since the 
cylindrical length is zero. Doubling the storage time from one to two years 
will almost triple the minimum storage mass. Figure 17 shows the effect of 
varying the length of the cylindrical cryostat section on minimum total system 
mass for a storage time requirement of three years. Again zero cylindrical 
length would correspond to a spherical container. The cylindrical lengths are 
in meters and do not describe the particular shapes of the containers. An option 
will be installed in the computer program which will allow tank length-to- 
diameter ratios to be parameterized. Figure 18 shows the influence which the 
percentage of tank fill has on the total system mass for a container with a three- 
year storage requirement, and a length of 3.05 m. For this solution, there is 
assumed to be no net heat input by an experiment load. As expected, the figure 
shows the total system mass increasing very rapidly for percent tank fill ap- 
proaching zero. Figure 19 shows the effect which the ratio of the experiment 
heat load and the MLI heat load has on the minimum total system mass, again 
for a storage time of three years and a cylindrical length of 3.05 m. The experi 
ment mass and volume are neglected. For Configuration 1, the experiment 
heat load and the penetration heat loads can be treated exactly the same. There- 
fore, since the experiment mass and volume are neglected, this figure is even 
more appropriate for penetration heat leaks since the mass of the penetrations 
are normally small compared to the total system mass. For the figure, 0., 

100., 200. and 300% penetration heat load corresponds to 0., 1, 2 and 3 times the 
MLI heat load, respectively. 

Configuration 2 is similar geometrically with Configuration 1 except that Con- 
figuration 2 has a vapor cooled shield (VCS) embedded within the MLI, as is 
shown in Fig. 20. The method of solution of this configuration is to initialize 
r j and the VCS temperature (T 2 ) and use the equation described for Configura- 
tion 1 to calculate r^ and the associated masses and heat fluxes. The rate of 
the cryogen vapor mass flow is calculated based on the initial mass of cryogen 
divided by the storage time, which assumes a steady state solution. The cool- 
ing capacity of the' VCS can then be calculated ~basedon“the'' shield" temperature"’ 
and the cryogen temperature as: 


q VCS ■ .^shield 


m 


C 

cr yo p 

vapor 


cryo 

vapor 


(T shield " T cryo ) 


where r) s hield transfer efficiency of the shield. The heat trans- 

ferred from the outer radius to the VCS, r 3 to r 2 * is then the heat transferred 
between the VCS and the cryogen tank, r 2 to rj f . plus the heat absorbed by the 
VCS, i.e., ’ ^ 


q 2-3 q l-2 +q VCS 
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Fig. 20 - Geometry for Configuration 2 


which must also be the same as the sum of the individual modes of heat transfer 
2 * 


from r ^ to r 9 , l .e . , 


q 2-3 ~ q MLI, , + q 

2-3 ^pene^ ^ 


where the calculation of q will be as described earlier. The heat 

P ene 2-3 


“transferred through “the MFI frorn r" to r^ is then 


1 MLI 2 _ 3 = q 1 - 2 + q VCS ' q pene 


Also, it is known that 


*MLI 


= 2jr 


2-3 


2-3 
T_ 


2 r 2 r 3 


(r 3 - r 2 ) + in(r,/r,) 


y z 


f 


k mli dT 


J T. 
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The prior two equations are then used to determine r%. When r3 is found for 
a particular VCS temperature, and after determining the mass of the vapor 
cooled shield from 


M VCS = Zp MAT r 2 + 2r r 2 L) thlck wall 

the mass of MLI between r 2 and r$ is calculated, and the total system mass 
is determined from 

, M, . . = + M wt t + 

„ total Config.l VCS MLI^ 3 pene^ 3 


which is stored. The VCS temperature is incremented and the above pro- 
cedure repeated until the minimum total system mass is determined. The 
above equations and procedure have been programmed into a computer routine 
.which is in the final stages of development. 

Summary of Assumptions Incorporated in Math Models : In, order to .construct 
a convenient set of parametric curves delineating the influence of configuration 
design on the mass of liquid helium systems, several simplifying assump- 
tions have been made which,' it is thought, do not severely compromise the 
validity of the mathematical models. These assumptions are summarized in 
« Table 4 for the basic design configurations. The assumptions made for Con- 
I figuration 1 are common.to all configurations. Configurations 2,3 and 4 re- 
quire additional assumptions. 

\ A preliminary matrix for the parametric study of liquid helium storage is „ 
.shown in^Table 5. ‘ The parameters considered to be most important are: 

^ . <■ . ' ' _ ‘ 

• Storage time r ’ ^ 

.. m Penetration conductance 

• Experimental heat load 

J o • Tank L/D 

* *\ ' 

These variables have been given emphasis in both range of parameter variation 
and in number of computer runs in the analytical investigation of their influ- 
ences on the minimum system mass required for liquid cryogen storage. The 
* parameters that will be varied to a lesser extent are: 

» “ ; t 

* ‘ " . • Tank ullage (includes experimental volume) 

,• Vapor shield efficiency (applicable to 
Configuration 2) • 

# Thermal short efficiency '(Configurations 2, 

3 and 4) 
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Table 4 

SUMMARY Of ASSUMPTIONS FOR PARAMETRIC STUDY 


Configuration 

Assumptions 

i 

One-dimensional heat transfer along 
radius 

• Cryostat wall is fully wetted 

• No localized boiling 

• No lateral temperature or 
property gradients in MLI 

Penetration conductance is constant 

Radial temperature gradient in tank 
wall is small 

Penetration channels are not vapor 
cooled > 

Weight of supports, vents, electrical 
leads, and experiment are not con- 
sidered in total system mass 

2,3.4 . 

Constant temperature shield, shroud, 
fin 

Penetration channels are cooled only 
at thermal shorts 

Penetration conductance is constant 
between external surface, shorts and 
c r yog en 


The following parameters will be fixed for analytical, study:. 

• External surface temperature of cryostat; 530°R 

• Liquid helium temperature: 5°R (saturated con- 
ditions) 

• MLI: superfloe at 11.8 layer/cm 

• Tank material: aluminum 
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Table 5 

PREUMINARY PARAMETRIC MATRIX FOR LIQUID CRYOGEN STORAGE STUDY 



Principal Matrix (192 rune) 


Delineate Influence of Ullage 
(12 rune) 


100 100 Principal Matrix (192 runt) 


75 100 Delineate Influence of Shield 

.. Efficiency (6 runt) 


Delineate Influence of Thermal 
Short Efficiency (9 rune) 


100 Principal Matrix (72 rune) 


Delineate Influence of Thermal 
Short Efficiency (9 rune) 


Principal Matrix (72 runs) 


Influence of Thermal Short 
Efficiency (9 runs) 


Principal Matrix (72 runt) 


Influence of Thermal Short 
Efficiency (9 runt) 


Principal Matrix (192 runt) 


Influence of Thermal Short 
Efficiency (9 runt) 
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Nomenclature 


TAP 

P 

P 

o 

T 

T 

o 

l/d 

Re 


pressure at time t 
initial pressure 

temperature or period of data record 
initial temperature 
length/diameter ratio of tube 
acoustic Reynolds number 


PSD 

x(t), y (t) 


x , 
i 


y 




x, y 



z 

1 




arbitrary functions of time 
discrete values of x(t), y(t) 
mean values 
mean square values 
standard deviations 
statistical variance 
complex function x^ + jy^ 

Fourier transform of z 
complex conjugate of ' 

power spectral density function. .. _ 
phase angle 
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Long Term Storage 


A 

area 

c 

penetration conductance 

c ; c 

P p 

^cryo *cry 

vapor 

specific heat of cryogen vapor at constant 
pressure 

k mli 

thermal conductivity of multilayer insulation 
perpendicular to the grain 

*K 

pene 

thermal conductivity of penetration channel 

L 

length of cylindrical section of the cryogenic 
storage tank 

L 

pene 

length. of .penetration channel 

^Config • 1 

mass of Configuration 1, if total system mass 
of a cylindrical tank with spherical end caps 
surrounded by multilayer insulation*, including 
cryogen, for specified conditions detailed in text 

M MLt 2 . 3 

mass of multilayer insulation between layers 
defined by radii r^ and r^ 

M 

P ene 2 . 3 

mass of penetrations through 

^total 

total system mass 

M vcs 

mass of the vapor cooled shield 

m , m 

cryo cry 

mass flow of cryogen vapor (boiloff) 

vapor 

__ — - ----- — ~ ^ ^ ~ ” 

q l-2 - - - 

heat transferred between radii r^ and r^ 

q 2- 3 

heat transferred between radii r^ and r^ 

q MLI 

heat transferred through multilayer insulation 

’ MLI 2.3 

heat transferred through multilayer insulation 
bounded by radii r^ and r^ 

qpene 2-3 

heat transferred through penetrations bounded 
by radii r^ and r^ 



Page 39 

LMSC-HREC PR D390395 
7 November 1974 


q l 


heat transfer rate from ambient environment to 
shield or shroud 


q 3 

q vc's ,q v.s. 


heat transfer rate to shield or shroud 

heat transfer rate from shield or shroud to 
c r yog en 

heat absorbed by the vapor cooled shield 
inner radius of cryogen tank 


radius for the vapor cooled shield 


outer radius for Configuration 2 (cryogen tank 
with a vapor cooled shield) 


amb 

T T T 
T 2’ 3 


cry 

thick 


wall 


ambient temperature 

temperature for r^,r ^ and r^ respectively 

temperature of stored cryogen 

thickness of walls such as the tank wall and 
the vapor cooled shield wall 


T 

e 

T sh 

^shield’ 17 V.S . 
Vs. 

^mat 

T P~ 


vapor temperature at exit of vapor cooled shield 

temperature of vapor cooled shield 

vapor cooled shield heat transfer efficiency 

thermal short efficiency 

material density 

temperature~at thermal short 
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New Technology 

The contract is being monitored continuously for new technology efforts that 
would be reportable under the contract clause entitled "New Technology." . 
Any pertinent items will be reported under that clause. 


Very truly yours, 

LOCKHEED MISSILES & SPACE COMPANY, Inc 

Juan K. Lovin, Supervisor 
Thermodynamics & Structures Section 
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